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ИССЛЕДОВАНИЕ ДИНАМИЧЕСКОЙ МОДЕЛИ 
АЭРОДРОМНОЙ УБОРОЧНОЙ МАШИНЫ 


Разработана математическая модель аэродромной уборочной машины, снабжённой 
снегоочистительным отвалом и вращающимся барабанным щёточным устройством. 
Произведен анализ этой модели и оптимизация её параметров с целью обеспечения 
качественной очистки аэродромного покрытия. 

Ключевые слова: динамическая модель, уборочная машина, поверхность аэро- 
дрома. 


Введение. Для обеспечения качественной очистки аэродромного покры- 
тия уборочная машина должна работать в режиме копирования поверхно- 
сти этого покрытия. Это гарантирует безотрывное движение рабочих орга- 
нов и выполнение ими заданных функций. Важным моментом построения 
математической модели машины является выбор эквивалентной схемы, 
отражающей степень идентичности модели реальному объекту. 
Постановка задачи. Задачей исследования является создание математи- 
ческой модели машины, эквивалентной реальному объекту, анализ этой 
модели и оптимизация её параметров с целью обеспечения качественной 
очистки аэродромного покрытия. 

Решение задачи. На рис.1 представлена модель аэродромной убороч- 
ной машины, которая состоит из тягача и шарнирно присоединённого к 
нему полуприцепа. Впереди машины установлен снегоочистительный отвал 
для предварительной очистки покрытия от снега. На полуприцепе машины 
установлен параллелограммный механизм с вращающимся барабанным 
щёточным устройством, производящим окончательную зачистку поверх- 
ности аэродрома от снега. 

Для обеспечения качественной очистки аэродромного покрытия 
уборочная машина должна работать в режиме копирования поверхности 
этого покрытия. Щёточное устройство (дальше - барабан) закреплено на 
заднем звене параллелограммного механизма, опирающегося на пневмати- 
ческие колёса, которые для обеспечения копирования должны соприка- 
саться с поверхностью аэродромного покрытия. 
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Рис. 1. Схема уборочной машины 


При составлении модели использованы уравнения Лагранжа с не- 
определёнными множителями с дополнительными голономными связями. 


а от ОГ 0П дФ хх ОЕ 
- +, -^=0О,, (1) 
4 04, 04, 04, 04, “—’ 904, 


где д - обобщенные координаты системы; Ги /7- соответственно кинети- 
1 





7 


ческая и потенциальная энергия системы; Ф - диссипативная функ- 
ция системы, О, — обобщенные силы системы; й — дополнительные 


голономные связи копирования, налагаемые на систему. 
В статье приняты следующие обозначения: [,— линейные размеры; 
®,— высоты неровностей поверхности; ги, го, гз, г. — соответственно коэф- 


фициенты сопротивлений шин ведущих и управляемых колес машины, си- 
стем навески и опорных колёс отвала и барабана; С;, С», Сз, С. — соответ- 
ственно коэффициенты жёсткости шин ведущих и управляемых колёс ма- 
шины, систем навески и опорных колёс отвала и барабана; Р., Е’, Е.- про- 


екции силы сопротивления Ё’, приложенной к отвалу на соответствующие 
оси координат; М. — момент сил сопротивления, приложенных к вращаю- 
щемуся барабану; М,, - вращающий момент, приложенный к барабану; 
т, — масса тягача; т, — масса полуприцепа; т. — масса отвала; т, — 
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масса рамы барабана; т, — масса барабана; Л„,Л,, — центральные мо- 
менты инерции тягача относительно соответствующих осей; Л, Л,, — цен- 
тральные моменты инерции полуприцепа относительно соответствующих 
осей; 7.„, /., - центральные моменты инерции отвала и рамы ба- 
рабана относительно оси ОХ; Л, — осевой момент инерции барабана отно- 


сительно главной центральной оси О.У; Л, — осевой момент инерции ба- 
рабана относительно главных центральных осей О.Х, и О.7,; й - смеще- 
ние центра инерции барабана от оси вращения; © - угол отклонения 
главной центральной оси от оси вращения барабана; у, и у, -— соответ- 


ственно углы установки отвала и барабана к поперечной оси машины. 
Движение колебательной системы рассматривается относительно 
инерциальной системы координат, движущейся поступательно, прямоли- 
нейно и равномерно, направление оси Ох которой совпадает с направле- 
нием вектора скорости движения машины. 
Введём обобщённые координаты 4,( }=1..9) следующим образом: 


4, = <, - вертикальное колебание тягача и полуприцепа; а, = ф, — угловые 
колебания тягача и полуприцепа относительно поперечной оси; 
4. =$, - угловые колебания тягача относительно продольной оси; 
4. =ф, - угловые колебания полуприцепа относительно продольной оси; 
4; ==, - вертикальные колебания отвала; а, =‹ф, - угловые колебания 
отвала в вертикальной плоскости; а, = <, — вертикальные колебания бара- 
бана; 4, =ф; - угловые колебания барабана в вертикальной плоскости; 
4 =ф - угол поворота барабана относительно собственной оси вращения. 


В режиме копирования рабочими органами машины (отвалом и ба- 
рабаном) на систему накладываются дополнительные голономные связи — 
связи копирования, которые вынуждают рабочие органы двигаться по по- 
верхности аэродромного покрытия без отрыва: 

= о + 2. ; аа —Ф 5. р 
Е = 8. = < + (ф. +; 4; Е = 8 =5-(ф.+Ф5Ь. 
Используя связи копирования, исключаем из системы четыре обоб- 


щённые координаты: 
_ 85186. 85 655 _ 871 88. 87 83 


› Ех: ‚ = › о абЗЕИЫ . 
5 2 21. $; 5 5 Ф; 21, Ф: 


Подставив выражения для кинетической, потенциальной энергий 


<> 


системы, диссипативной функции и обобщённых сил в уравнения Лагранжа 
(1), после преобразования получим систему дифференциальных уравне- 
ний, адекватно описывающую динамику рассматриваемой системы: 


1-1 
(т, +т, ) 5 +2(*, т а +2(С, и С.) = + 2(151 — 1/5) ф, +2(С>14 - С), р 
п (2, + $2) — 1 (83 +84)- С1(в + 82) С5 (83 +84)=0. 
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1Е2 
ве. +.,, +1.) +2(1 +516, +2(С/ +С р, 3 ф› + 2(51, ть) я 
+2(Сы, - С) = +, ет 5 + 71, (8, +8, )- "51, (8, +8.) + С. (в, +8,)- 


— Си (в, +8.) =0. ь 
7=3 





й р 1 | : 
‚$ +215 Ф, +2С6, с» ыы -&)+ 


10 
2 


+71 (85 -&;)+® 1—8, - &)- С. (в, -&)-Аь + 1.10 Е =0. 
10 


1-4 
1$, + (1+ 1,6 бсо$у, )со$(ф)- 
а. го 83 соз? ( ‚)5°ф соз(ф)зт (=)- (7, + 1.6 ф’ соз( ‚5 (Ф)+ 


8 





м 8—8 8 соз? ( ‚6? соз? ()+ (1. соз( ‚+. + 
51. 57. 


+ 5 (8 -&,)+ С (8. -&,)- №, + 7, = 0. 


а 





т, 85—86 ДЛЕ =0, 
2 
7=6 


:. .. 85-8 . Е 
Уз. Фу — Уз фо + Ух > о [в 851 - `, 
10 10 











85-8 
В 851 - ` ЯщЦо + Я>о — Е =0. 
10 
ИЕР 


т;й с0(ф) ф— туй эт (ф) 6? + (т, + т, - 21. (8. +8.) 

-2С. (в. +8:)-№ 1, =0. 

7=8 

(7 + 75) бсоз(у, )соз(ф)в 2 ВТ 8 сз? (9,5 ?ф зи(ф)соз(ф)+ 
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183 
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+. —_ 8 соз? (у,)- (7, +; )бсоз(у 6 зт(ф)+ 
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т 87 83 
21 
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758 соз" (у, )6° соз? (ф) + Ла, “:- ЛЬ ЕГО 
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1=9 
(т.т? +7 + 1.5?) ф+т.йсоз(ф) 8 св (7, }- 8 6 соз(у, )соз(ф) + 


8 





8 

Для проведения исследований использована система МАТЕАВУ.О.1 
и её подсистема моделирования динамических процессов 5итштпк. Реше- 
ние системы уравнений 1, 2, 3, 5, 6 производим с использованием блока 
ОШегепца! ЕдчаЧоп ЕЧКог. Схема модели приведена на рис. 2, а результаты 
решения на рис. 3. 
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Рис. 2. Схема модели 


На рис.3‚а кривая отображает вертикальные колебания машины; 
на рис.3,б,в кривые — угловые колебания машины относительно продоль- 


ной и поперечных осей; на рис.3‚г,д кривые отображают значения ФД, и 4, 


- реакций опорной поверхности на копирующие колёса отвала. 

Результаты исследования показывают, что неровности копируемой 
поверхности вызывают колебания рабочего органа - отвала. Эти колебания 
приводят к изменению величины опорных реакций на копирующих колёсах 
отвала. Минимальное значение опорной реакции значительно меньше ста- 
тической нагрузки(от силы тяжести) на опорных колёсах. Это гарантирует 
безотрывное движение рабочего органа и выполнение им заданных функ- 
ций. 
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Рис.3. Результаты решения системы уравнений 


Для оптимизации параметров щёточного барабана использованы 
уравнения 7, 8, 9. В результате решения уравнения 9 получена зависи- 
мость угловой скорости & вращения щёточного барабана от времени. Вид 
этой зависимости показан на рис. 4. Прямая ОА (см. рис.4) соответствует 
разгону ротора без рабочей нагрузки, точка А соответствует моменту 
включения рабочей нагрузки (щёточный барабан опускается в рабочее по- 
ложение), кривая АВ характеризует переходной процесс, а далее вращение 
барабана происходит с установившейся угловой скоростью. 

Привод щёточного барабана осуществляется гидромотором, имею- 
щим крутую характеристику (большим изменениям нагрузки соответствуют 
малые изменения угловой скорости). Значение вращающего момента взято 
из характеристики гидромотора. Значение момента М. сопротивления 
зависит от величины угловой скорости @ и имеет вид: 

М 
где М - постоянная составляющая момента сопротивления; к- постоян- 


ный коэффициент. По результатам исследований М=5 Нм, 
к = 17 Нмс. 


се ЕМ+ ко, 
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0 1 2 3 4 5 Б 7 8 3 10 #5 
Рис.4. Зависимость угловой скорости @ вращения от времени # 


Из уравнений 7 и 8 находим значения д, и Д,, характеризующие 


величину динамических реакций копирующих колёс щёточного устройства. 
Результаты решения показывают, что величина опорных реакций копиру- 
ющих колес определяется параметрами, характеризующими степень ба- 
лансировки барабана: К - величина смещения центра масс барабана от оси 
его вращения, 6 - угол отклонения главной центральной оси барабана от 
оси вращения. 

На рис. 5,а приведены усилия на копирующих колёсах при сбалан- 
сированном барабане (п=0, б =0) для разных скоростей движения маши- 
ны. Величина динамических реакций, вызванных неровностями копируемой 
поверхности, возрастает с увеличением скорости движения машины. 
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Рис.5. Динамические реакции копирующих колёс 


Размах колебаний динамических реакций при максимальном значе- 
нии амплитуды возмущающего воздействия составляет 4000Н, что значи- 
тельно меньше статического усилия - силы тяжести, приходящейся на 
опорное колесо. Однако при недостаточно сбалансированном барабане эти 
усилия могут достигать больших значений, превышающих величину силы 
тяжести. Необходимым условием копирования является положительное 
значение реакции на копирующем колесе. Исходя из этого условия могут 
быть определены максимально допустимые значения параметров ВН и дб. 

На рис.5,б приведено значение усилий на опорных колёсах не пол- 

ностью сбалансированного барабана: 6 =0.4 градуса, Н= 0.004м. Эти зна- 
чения вычислены для скорости движения машины 40 км/ч. При таких зна- 
чениях дисбаланса барабана, как видно из графика (см. рис.5,6), условия 
копирования выполняются (минимальное значение опорной реакции поло- 
жительное). При превышении указанных пределов одного из параметров 
били В минимальное значение опорной реакции становится меньшем ну- 
ля и условия копирования нарушаются. 
Выводы. Получены уравнения, являющиеся математическим описанием 
динамики рассматриваемой колебательной системы. На вход системы по- 
даются возмущающие воздействия со стороны неровностей очищаемой по- 
верхности и сил сопротивления, а также воздействия от неуравновешенно- 
сти барабана. На выходе присутствуют неопределённые множители Ла- 
гранжа, значения которых использованы для математического формули- 
рования критерия оптимизации системы. 

Проведен анализ математической модели машины и оптимизация 
её параметров. В результате анализа этой модели установлено, что снего- 
очистительный отвал выполняет заданный технологический процесс очист- 
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ки на скорости движения машины до 40 км/ч. При этом величина динами- 
ческих реакций на копирующих колесах, возрастающая с увеличением ско- 
рости движения машины, остаётся в допустимых пределах. Выполненная 
оптимизация параметров исследуемой модели позволила установить пре- 
дельные значения параметров дисбаланса щеточного барабана для работы 
на скорости до 40 км/ч. 
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